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Abstract 

Based on the 7-ray data of blazars in the third EGRET catalog and radio data at 5 GHz, we studied 
the correlation between the radio and 7-ray luminosities using two statistical methods. The first method 
was the partial correlation analysis method, which indicates that there exist correlations between the 
radio and 7-ray luminosities in both high and low states as well as in the average case. The second 
method involved a comparison of expected 7-ray luminosity distribution with the observed data using the 
Kolmogorov-Smirnov (KS) test. In the second method, we assumed that there is a correlation between the 
radio and 7-ray luminosities and that the 7-ray luminosity function is proportional to the radio luminosity 
function. The KS test indicates that the expected gamma-ray luminosity distributions are consistent with 
the observed data in a reasonable parameter range. Finally, we used different 7-ray luminosity functions 
to estimate the possible 'observed' 7-ray luminosity distributions by GLAST. 

Key vi^ords: galaxies: active — quasars: general — gamma rays 



1. Introduction 

Although the correlation between the radio and 7-ray 
luminosities of blazars has been widely studied, two op- 
posite conclusions have been made. Based on the data 
observed by EGRET at different observation times, a 
moderate correlation between the radio and 7-ray lu- 
minosities has been obtained (e.g. Stecker et al. 1993; 
Padovani et al. 1993; Salamon, Stecker 1994; Dondi, 
Ghisellini, 1995; Stecker, Salmon , 1996; Miicke et al. 
1997; Mattox et al. 1997; Zhou ct al. 1997; Fan et al. 
1998; Cheng et al. 2000). However, Miicke et al. (1997) 
performed Monte Carlo simulations of selected data sets 
and applied different correlation techniques in view of a 
truncation bias caused by sensitivity limits of the surveys. 
In their analysis, they used radio and 7-ray luminosity 
functions given by Dunlop and Peacock (1990) and Chi- 
ang et al. (1995). After a careful analysis, they claimed 
that there is no correlation between the radio and 7-ray 
luminosities for simultaneously observed radio and 7-ray 
data. 

Many models have been proposed to explain the origin 
of the blazar 7-ray emission, including synchrotron self- 
Compton (e.g. Maraschi et al 1992), inverse Compton 
scattering on photons produced by an accretion disk (e.g. 
Dermer et al. 1992; Zhang, Cheng 1997), or scattered by 
ambient material, or reprocessed by broad-line clouds ( 
e.g. Sikora et al. 1994; Blandford, Levinson 1995; Xie et 
al. 1997), synchrotron emission by ultra-relativistic elec- 



trons and positrons (e.g, Ghisellini et al. 1993; Cheng 
et al. 1993), and electromagnetic cascade by collision 
of ultra-relativistic nucleons (e.g, Mannheim, Biermann 
1992; Mannheim 1993; Cheng, Ding 1994). However, 
there is no consensus yet on the dominant emission pro- 
cess. 

Correlation analye between different wavebands are 
important for us to understand the 7-ray emission from 
blazars. For example, if radio emission is produced by 
high-energy electrons in the jet via synchrotron radia- 
tion, it is likely that these electrons contribute part of 
the 7-rays due to inverse Compton scattering. In this 
sense, it seems reasonable to expect some correlation be- 
tween the radio and 7-ray luminosities. Here, we revisit 
the correlation between the radio and 7-ray luminosities 
of blazars using the observed 7-ray data in the third cat- 
alog of EGRET (Hartman et al. 1999) and radio data at 
5 GHz (see table 1), since the data in these two bands 
are presently complete. In section 2, we list the observed 
data and discuss a correlation analysis based on the use 
of a partial correlation analysis. Further, we estimate 7- 
ray luminosity distribution of the observed blazars at 100 
MeV and compared it with the observed data in section 
3. Finally, we give our conclusions and a brief discussion 
in section 4. 

Throughout the paper we use a Hubble constant, Hq, 
of 75 km s~^ Mpc~^ and the deceleration parameter, go, 
of 0.5. 
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Fig. 1.. Correlations of the 7-ray luminosity with radio luminosity. 
(A) low state. (B) high state. (C) average case. The differential 
luminosity is used in a specific frequencies u, where i/ is 5 GHz 
for the radio and 100 MoV for the 7-rays. 



2. Observed Data and Correlation Analysis 

The thhd EGRET catalog includes 66 positive and 25 
marginal detections of 7-ray loud blazars above 100 MeV 
(Hartman et al. 1999). We list 67 of the identified blazars 
(including Markarian 501) and 25 marginal blazars in 
table 1. For each source, its 7-ray fluxes above 100 MeV 
in both the high and low states are the maximum and 
minimum observed data in the third EGRET catalog. 
The average 7-ray flux is the value labeled by P1234 in 
the third catalog. For radio data at 5 GHz, part of our 
sample were taken from those in both high and low states 
compiled by Cheng et al. (2000). The average radio data 
at 5 GHz were taken from a database provided by the 
University of Michigan Radio Astronomy Observatory. 



From table 1, we can obtain that the mean value of 
the redshifts in our sample is 0.96 ± 0.06 and the average 
value of integrated spectral index of 7-rays is 1.22 ± 0.01 
using a weighted average method. In order to analyze 
the correlation between the radio and 7-ray luminosities, 
we made a K-correction for the observed fluxes listed in 
table 1 by using 



(1) 



where z is the redshift and a is the spectral index at 
the frequency v. We used a = in the radio band and 
a — 1.22 in the 7 energy range when the spectral indices 
of the sources were not known. When the redshift was 
not known, we used the mean value of our sample, i.e. 
z = 0.96. We now consider the differential luminosity, 
which is given by 



Eif 



(2) 



where ai is the spectral index, i = R is for radio and 
i = 7 for 7-rays. We chose E^i = 100 MeV. If the flux 
density at a certain energy and the spectral index are 
known, the differential luminosity per unit solid angle is 
given by P,{ = S,{Ei{)[Rlr'^{l + z)^], where S, is the 
radio flux density for z = R or 7-ray flux density for 
i = 7, z is the redshift, r is the co- moving coordinate of 
the source, Rq is the present cosmological scale factor, 
and Ror = {2c/Ho)[l -{1 + z)-'^/% Therefore, for radio 
data. 



Pi 



Rf 



(1 



10^'" (t) [l^(l + ^)-^/¥ (3) 
2)""+^^Robs(PRf) WHz-1 sr-1 , 



where Sji obs(^'Rf) is in units of Jy. For 7-ray data, 

P,f « 6.6 X 10-« )' [1 - (1 + zr'/^a, (4) 

(1 -t- z)"-'+iFobs(> 100 MeV) W Hz"! sr'^ , 

where -Fobs (> lOOMeV) is the integrated 7-ray flux above 
100 MeV in units of lO"'^ cm'^ s'^ observed by EGRET 
and -|- 1 is differential spectral index. 

Using the observed data in table 1 and equationsq (jj) 
and (^, we obtained the differential luminosities at 5 
GHz and 100 MeV. We could thus consider the correla- 
tion between the radio and 7-ray luminosities in both the 
high and low states and for the average case. In the low 
state, the correlation between luminosities gives 

logioPRf = (13.7 ± 0.7) + (0.75 ± 0.04)logioP7f -(5) 
In the high sate, we have 

logioPRf = (15.0 ± 0.7) + (0.67 ± 0.04)logioP^f .(6) 
For the average case, we have 

logioPRf = (14.1 ± 0.7) + (0.74 ± 0.04)logioP^f .(7) 



The Radio and Gamma- Ray Luminosities of Blazars 
Table 1. Blazar sample*. 
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Table 1. Continued 



Source*^ 


z 




s 


H(Jy) 




s 


R (Jy) 




S'S'(Jv)'' 


r J 




Fi 




r J 






7 




0459-1-060* 


1, 


.106 




















3.40± 1.82 


1.19± 


0.32 


0.61+ 0.21 


1 


.36+0, 


.23 


0616-116* 
























4.79± 1.63 


1.49± 


0.82 


1.09± 0.34 


1 


.67+0, 


.27 


0738+5451* 





.723 




















4.21± 0.83 


1.14± 


0.63 


1.11+ 0.24 


1 


.03+0 


.10 


0850-1213* 


0, 


.566 




















4.44+ 1.16 


1.40+ 


0.44 









.58+0 


.21 


1243-072* 


1 


.286 










1 


1 A .,1, 


U, 


.04 




4.41± 2.96 


0.60± 


0.25 


0.98± 0.21 


1 


.73+0, 


.11 


1334-127* 


0, 


.539 




















2.02± 1.16 


1.14± 


0.38 


0.55+ 0.19 


1 


.62+0, 


.26 


1725-1-044* 


0, 


.296 










1 


f>rr 1 


U, 


.04 




3.02± 1.88 


1.33± 


0.61 


1.79+ 0.41 


1 


.67+0, 


.16 


1759-396* 






















2.29± 0.06 


14.57+ 4.89 


1.75+ 


0.47 


0.98+ 0.29 


2, 


.10+0, 


.24 


1908-201* 






2 




U. 


,03 


2 


1 n 1 


U, 


,05 


0.65± 0.03 


o. / lit 


Z.UO 


1.49+ 


0.34 


1.75+0.27 


1 


.39+0, 


.10 


1936-155* 


1, 


.657 


1 




U. 


.08 







U, 


,03 






l.oD 


0.76± 


0.30 


0.74 




2, 


.45+0, 


.90 


2022-077* 


1 


.388 




















7 4^+ 


1.04 


2.18± 


0.38 


2.12+ 


0.35 


1 


.38+0, 


.10 


2320-035* 


1 


.411 




















Q QO 1 


i.Ul 


0.82+ 


0.44 


0.60 










2351-456* 


1 


.992 




















4.28+ 


2.03 


1.18+ 


0.52 


1.43+ 


0.37 


1 


.38+0, 


.22 


0234+285 


1 


.213 


4 


QV 1 


u. 


,06 


1 


/I 1 


U, 


,04 


2.64± 0.07 


3.41+ 


1.16 


1.09+ 


0.44 


1.38+ 


0.26 


1 


.53+0, 


.13 


0506-612 


1 


.093 


1 


rrf\ 1 


U. 


.08 


1 


rrt 1 


fl 

u. 


,08 




2.88± 


1.15 


0.64± 


0.12 


0.72+ 


0.17 


1 


.40+0, 


.15 


0521-365 





.055 


8 




u. 


,17 


6 


.oz± 


U, 


,16 


8.06± 0.04 


3.19± 


0.72 


1.95± 


0.44 


1.58+ 


0.35 


1 


.63+0, 


.26 


0804+499 


1 


.433 


1 


CkA 1 


U. 


,13 





004- 


fl 

u. 


,04 




1.51± 


0.61 


0.83± 


0.39 


1.07+ 


0.25 


1 


.15+0, 


.24 


0917+449 


2 


.180 


2 


O A 1 


U. 


,06 


1 


nn 1 


fl 

u. 


,02 




3.35± 


1.30 


1.14+ 


0.33 


1.38+ 


0.20 


1 


.19+0, 


.08 


1127-145 


1 


.187 


5 


/JO 1 


U. 


,17 


2 


Qtr 1 


fl 

u. 


,10 


3.63± 0.02 


6.18+ 


1.80 


1.08+ 


0.59 


0.99+ 


0.24 


1 


.70+0, 


.20 


1313-333 


1, 


.210 


1 


A{±. 


U. 


,03 


1 


.zi± 


fl 

u. 


,18 




3.18± 


1.90 


1.62± 


0.53 


1.46+ 


0.25 


1 


.28+0, 


.11 


m 1 0-1-041 * 


0, 


.637 


2, 


.07± 


0. 


,09 


1 


.22± 


0, 


,06 




2.03± 


0.58 


1.34+ 


0.43 


0.51+ 


0.27 


1 


.63+0, 


.41 


0130-171* 


1, 


.022 










1 


.OOzb 


0, 


,07 




1.38+ 


0.68 


0.92+ 


0.52 


1.16+ 


0.30 


1 


.50+0, 


.17 


0415+379* 


0, 


.049 


8, 


m± 


0. 


,19 


5 


.81± 


0, 


,04 


6.44± 0.09 


6.02± 


1.71 


1.02+ 


0.31 


1.28+ 


0.26 


1, 


.59+0, 


.20 


0537-286* 


3, 


.110 


1 


.19± 


0. 


,10 


1, 


.02± 


0, 


,07 


1.16± 0.08 


3.50± 


1.18 


0.96+ 


0.49 


0.69+ 


0.27 


1 


.47+0, 


.36 


0539-057* 


0, 


.839 










1, 


.56ib 


0, 


,08 








6.65± 


1.95 


1.00 










0803-^5126* 


1, 


.140 




















2.34± 


1.37 


0.99+ 


0.26 


0.87+ 


0.24 


1, 


.76+0, 


.22 


1011+496* 


0, 


.200 




















0.80+ 


0.30 


0.44+ 


0.24 


0.48+ 


0.14 





.90+0, 


.18 


1055+567* 


0, 


.410 




















1.61± 


1.01 


0.65± 


0.16 


0.50+ 


0.14 


1, 


.51+0, 


.28 


1237+0459* 






















1.52± 0.86 


0.62± 


0.21 


0.65± 


0.15 


1.48+0.27 








1324+224* 


1, 


.400 


1 


.66± 


0. 


,17 


1 


.29± 


0, 


,15 


1.52± 0.07 


6.84± 


2.26 


0.95± 


0.27 


0.52+ 


0.16 


1 


.58+0, 


.16 


1504-166* 


0, 


.876 


2 


.72± 


0. 


,07 


2 


.12± 


0, 


,08 


2.39± 0.03 


3.32± 


1.03 


1.65+ 


0.63 


0.88 










1514-241* 


0, 


.042 


3 


.00± 


0. 


,05 


1 


.35ib 


0, 


,04 


2.23± 0.08 


3.72+ 


1.83 


1.85+ 


0.61 


0.84+ 


0.28 


1 


.66+0, 


.27 


1716-771* 
























4.05± 


2.31 


1.54+ 


0.57 


0.84+ 


0.40 


1, 


.74+0, 


.24 


1808-5011* 
























6.21± 


1.97 


0.73+ 


0.33 


0.59+ 


0.27 


1 


.93+0, 


.28 


2105+598* 
























3.32+ 


0.96 


1.62+ 


0.76 


1.98+ 


0.41 


1 


.21+0, 


.14 


2206+650* 
























3.08+ 


1.32 


1.86+ 


0.82 


2.44+ 


0.55 


1 


.29+0, 


.15 


2250+1926* 
























6.22+ 


2.15 


1.06+ 


0.38 


0.58+ 


0.28 


1 


.36+0, 


.35 


2346+385* 


1, 


.032 




















3.75± 


1.03 


0.85± 


0.36 


0.61+ 


0.32 


1, 


.47+0, 


.40 



: 7-ray flux is in units ofl0~^ cm~^ s~^ 

a: the data of source names which arc unmarked by stars arc taken from those compiled by Cheng et al. (1999), the radio 

data of the source name which are marked by stars in high state or low state are taken from Kiihr et al. (1981). 

b: the averaged radio data are taken from the database provided by the University of Michigan Radio Astronomy Observatory. 
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As pointed out by Padovani (1992), because our sam- 
ple was flux-limited, the luminosity was strongly corre- 
lated with redshift, would result in a spurious correlation. 
Therefore, we used a partial correlation analysis in our 
analysis, i.e. we examined the correlations between the 
radio and 7-ray luminosities excluding the dependence 
on the redshift. The correlation coefficients used for the 
analysis were the Spearman Rank-Order correlation coef- 
ficients (Press et al. 1992). In the low state (72 sources), 
the partial correlation coefficient is r w 0.26 (with chance 
probability P « 2.6 x 10~^), which was derived by using 
the fact that the correlation coefficients between log P~^[ 
and log Pjif, log Pjf and log z, and log Prif and log z are 
0.904, 0.972, and 0.903 respectively. In the high state 
(62 sources), the correlation coefficients between log Pj{ 
and log Prif , log Pjf and log z, and log Prif and log z are 
0.919, 0.967, and 0.919 respectively; the partial correla- 
tion coefficient is thus r w 0.30 (with chance probability 
P « 1.8 X 10^^). Therefore, we may conclude that there 
exist correlations between the radio and 7-ray luminosi- 
ties in both the high and low states. For the average 
case, from table 1, there are 51 sources which have aver- 
age observed radio data. The average 7-ray flux for each 
source is the value (labeled P1234) in the third EGRET 
catalog (Hartman et al. 1999). The partial correlation 
analysis gives r « 0.28 with a chance probability, P, of 
w 3.9 X 10^^, which is at least a marginal correlation. 
Figure ^ shows the correlation between the luminosities 
for the above three cases. 

3. Gamma-Ray Luminosity Distribution 

We now consider the 7-ray luminosity distribution in 
order to further discuss whether there is a correlation 
between the radio and 7-ray luminosities. Using a pre- 
existing radio luminosity function and assuming that the 
radio and 7-ray luminosity are related [the correlation 
coefficients are given by equations (^, (||) and (|^)], we 
could calculate the 7-ray luminosity distributions in the 
low, high states, and average cases following the method 
given by Stecker et al. (1993). As a comparison, when we 
calculated the 7-ray luminosity distribution in the aver- 
age case, we used 7-ray luminosity function given by Chi- 
ang and Mukherjee (1998). These expected distributions 
were compared to the observed distribution. In figure 
the normalized cumulative observed 7-ray luminosity 
distributions in the low, high states, and the average case 
are shown (dotted curves). 

If the 7-ray luminosity function (pj) is given, the total 
number of sources of the luminosity, P^f, seen at Earth 
is (Salamon, Stecker 1994) 

nobsd(logioP7f ) = /J^"'"^" 4^p3r2 (8) 

p^{P^i,z)dr d{\ogj^QPy[) , 




Fig. 2.. The Comparison of observed and expected cumulative 
distributions of the blazar 7-ray luminosities at 100 MeV. (A) 
low state. (B) high state. (C) average case. For the average 
case, the 7jLLay luminosity distributions calculated by using Eqs. 
(12) and (hj) are represented by solid, dashed (71 = 1.2 and 
72 = 2.2) and dot-dashed (71 = 2.9 and 72 = 2.6) curves re- 
spectively. Dotted curves represent the observed distributions. 
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Fig. 3.. The comparison of the differential 7-ray luminosity di stri - 
butions calculated by using Eqs. ( |l2| ) (dotted curve) and (|l4) 
(dashed curve) with the observed data (histogram). 
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where P^f is given by equation (|^) for a source with spec- 
tral index a^, and -RoJ^max is the maximum detectable 
distance which is found by using equation (|^) and sub- 
stituting the flux sensitivity, P^um, for Pobs- We can 
convert i?o^max to Zmax using the relation of Ror and z. 
In other words, the maximum redshift can be estimated 

by 



f/[(2c/i7o)'5^ lim] . 



(9) 



It should be noted that the co- moving density, p-y, is in 
units of Mpc~^x(unit interval of logigP-y)"^. 

Salamon and Stecker (1994) assumed that (i) there is 
a linear relation between the 7-ray luminosity (Pji) and 
the radio luminosity, i.e. P-yf = lO^Ppif and (ii) Pj = 
rjpB,, where 77 is a normalization factor and pr is a radio 
luminosity function. They found ^ = —9.8 by fitting 38 
EGRET sources. Here, we assume that Pnf = lO^P^f, 
where the values of a and b in the high, low states, and 
for average case were estimated by equations (||), and 
(0) respectively. The radio luminosity function is given 
by Dunlop and Peacock (1990) 



Pr(-Pr, z) 



1 


\ 1 


0.83 

+ 




1.96 1 




[Pc{z)\ 









,(10) 



where LogP^iz) = 25.26 + l.lSz - 0.28^^ Pr is differ- 
ential radio luminosity in units of W Hz~^ sr~^, and pR 
is co-moving density in units of Mpc^'^ x (unit interval of 
logioPR)"^. Using above assumptions and equation (p^, 
we have 

nobsd(logioP7t) = 2^77 (1)' /o'""" (11) 

lii±f^^PR(10'^P^^,z)dzd(logioP7f) , 

where Zmax is determined by equation (^^. Generally, 
once the Zmax is given, we can obtain the 7-ray luminos- 
ity distribution. Because we assumed a maximum cutoff 
redshift of z^ax — 5, the upper limit of equation ( p^ ) is 
min (zinax, 5)- In order to estimate the 7-ray luminosity 
distribution using equation (p^), we needed to know the 
7-ray flux sensitivity and to then determine two parame- 
ters (a, b). We adjusted a and b within the allowed range, 
which is limited by equation (||) in the low state or by 
equation (||)in the high state. We used a = 13.7, b = 0.76, 
ii„i(> 100 MeV)= 0.5 x IQ-'^cm-^ s~i in the low 
state, a = 14.0, b = 0.68, um{> 100 MeV)= 1.5x10^^ 
cm^^ s^^ in the high state, and a — 14.1, b — 0.74 and 
P7 iim(> 100 MeV)= 0.5 X 10"^ cm'^ s'^ for the average 
case. We used one sample KS test to compare the cumu- 
lative distributions with the observed ones. The max- 
imum deviation are 0.09 in the low state (92 sources), 
0.1 in the high state (91 sources) and 0.07 for the av- 
erage case (92 sources). Therefore, the null hypothesis 



■F 0.6 




13 14 15 16 17 
log (P )(WHz"'sr"') 



Fig. 4.. The expected differential distributions of average 7-ray lu- 
minosities of the blazars at 100 MeV. The threshold of GLAST 
is used (F^ iim{> 100 MeV) = 3 X lO'^ cm-^ s-lL_The 7-ray 
luminosity distributions calculated by using Eq. ma) and Eq. 
( p^ are represented by solid and dashed curves respectively. 



cannot be rejected at the > 80% confidence level. Com- 
parisons of the observed data and our results in the low, 
high states, and the average case are shown in panels (A), 
(B), and (C) of figure ^ respectively. 

For a comparison, in the average case, we could also 
estimate the average 7-ray luminosity distribution us- 
ing the 7-ray luminosity function of the EGRET blazars. 
The 7-ray luminosity function can be expressed as (Chi- 
ang and Mukherjee 1998) 



P7 EGRET 



(*) 



-72 



Pb 



e(P70 - Pb) 
Q{Pb 



(12) 



p 



70; 



where 71 « 1.2, 72 = 2.2±0.1, logPe « 14.9 W Hz"! sr"! 
(corresponding to Lb ~ 1-1 x lO''^ erg s~^) and 8(x) is 
a step function. P^q = P^/ f{z), where f{z) = (1 + z)^ 
with f3 — 2.7 ± 0.7. Using equation (^), we can rewrite 
equation 



as 

nobsd(logioP7f) = 47rlogJ0 (|^)^ /q' 



(13) 



(l-|-z)5/2 



-P7 EGRETP-yfdz (i(logi(,P^f ) 



where the upper limit of the integral is min(2:fnax, 5). We 
calculated the 7-ray luminosity distribution using equa- 
tion (|l^. In our calculation, we used the minimum value 
of the average 7-ray flux in our sample as the flux sensi- 
tivity, which is consistent with that used by Chiang and 
Mukherjee (1998), and (3 = 2.2. The maximum devia- 
tion is 0.12 for the expected 7-ray luminosity distribution 
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calculated by equation (p_4[). Therefore, the KS test indi- 
cates that the null hypothesis cannot be rejected at the 
> 90% confidence level. A comparison of the observed 
and expected cumulative 7-ray luminosity distributions 
is shown in the panel (C) of figure 

4. Conclusion and Discussion 

Using 7-ray data of the blazars in the third EGRET 
catalog (Hartman et al. 1999) and radio data at 5 GHz, 
we used two methods to consider the correlation between 
the radio and 7-ray luminosities in the high and low 
states as well as in the average case. We first made a 
partial correlation analysis for the radio and 7-ray lu- 
minosities. Our results indicate that there are, at least 
marginal, correlations between the radio and 7-ray lumi- 
nosities for the three cases. Second, we calculated the 7- 
ray luminosity distributions of the blazars by assuming a 
correlation between the radio and 7-ray luminosities (the 
correlation coefficients were determined by the observed 
data, see equations ( (||) and (|^) for both the low 
and high states as well as the average case) and using 
the 7-ray luminosity functions given by equation (^2|). 
We found that the expected 7-ray luminosity distribu- 
tions are consistent with the observed data by using KS 
test (see figure [2]) . Therefore, after using two statistical 
methods, we conclude that there is at least some corre- 
lation between the radio and 7-ray luminosities within a 
reasonable parameter range. 

We also used the 7-ray luminosity functions given by 
Chiang and Mukherjee (1998) to estimate the 7-ray lu- 
minosity distribution. It should be noted that in their 
calculation of the 7-ray luminosity function, Chiang and 
Mukherjee (1998) used the 1 Jy radio catalog of Kiihr et 
al. (1981) as an additional criterion to reduce the selec- 
tion effects resulting from an incomplete redshift sample 
of AGNs. Therefore, there are 34 AGNs in their sample. 
The excluded blazars either fail to meet the significance 
Hmit, (TS")^/^ = 4, or the 1 Jy radio flux hmit. In the cal- 
culation, they introduced a radio completeness function 
[see equation (7) of their paper] to reflect that the sam- 
ple is restricted to being in the Kiihr catalog. Therefore, 
although they did not explicitly assume a correlation be- 
tween the radio and 7-ray luminosities, they did include 
an assumption that all blazars must be observed at the 
radio band and that their radio energy fluxes must be 
greater than 1 Jy. Obviously, this assumption is a key 
factor by which the expected gamma-ray luminosity dis- 
tribution can roughly explain the observed data. In or- 
der to account for it, we note that Chiang et al. (1995) 
used the first EGRET catalog to obtain a gamma-ray 
luminosity function by modeling the optical evolution of 
their sources simultaneously with the gamma-ray evolu- 
tion. The obtained 7-ray luminosity function is the same 
as equation (O) with different indexes ( 71 = 2.9 and 



72 = 2.6). When we used the result of Chiang et al. 
(1995) to estimate the gamma-ray luminosity distribu- 
tion, the result was not consistent with the new EGRET 
catalogat all (see dot-dashed curve in the panel (C) of 
figure!]). 

From the KS test, we can see that the 7-ray luminosity 
distributions calculated by using equation (|l^) is better 
than that calculated by using equation (|lj). The former 
indicates a peak at logP-yf ~ 16.75, which is consistent 
with the observed data. The later indicates a wider peak 
from logPjf 15 to logP^f 17.5. In figure |3], we show 
a comparison of the expected differential 7-ray luminos- 
ity distributions with the observed data. Further, we use 
the possible threshold of GLAST, F-y - 3 x 10"^ cm'^ 
(Kamae et al. 1999), to expect the 7-ray luminosity 
distributions using equations (03) and (14), respectively. 
The results are shown in figure|4.|. Here, we have nor- 
malized the distributions to make the areas covered by 
the curves be the same. It can be seen that they have 
roughly the same shapes, but the former peaks are at 
logPyt = 15.6 and the later at logP.yf = 16.3. We expect 
that this difference will be examined by GLAST. 

From an observational point of view, there are some 
differences between BL Lacertae objects and flat spec- 
tral radio quasars (FSRQs), with the former showing 
no/weak emission lines, steeper X-ray spectra, and the 
later showing strong emission lines and flatter X-ray spec- 
tra. They are discussed separately in some cases (see 
Fossati et al. 1998; Ghisellini et al. 1998). But they 
both show a similar fiat radio spectrum and a similar 
GeV 7-ray spectrum. Since we only considered the radio 
and the 7-ray regions, we discussed them both together. 

Based on a statistical analysis, we have discussed the 
correlation between the 7-ray and the radio regions using 
the higher frequency radio data, the higher and lower 
states 7-ray and the radio regions (Fan et al. 1998; Cheng 
et al. 2000). We found that there is a closer correlation 
between the 7-ray emission and the high frequency (1.3 
mm, 230 GHz) radio emission for the maximum data 
than between the 7-ray and the lower frequency (5GHz) 
radio emission (Fan et al. 1998). When we revisited 
the multiwavelength correlation using the compiled data, 
we found that there are very strong correlations between 
Fx and Fq, and Fq and Fk in both the low and high 
states. However, a strong correlation between Fx and 
Fk exists only in the low state. There are also hints of 
an anti-correlation between F^ and Fx as well as Frmj 
and Fq, and a positive correlation between Fj and F/j. 
But the correlation between the 7-ray and the radio is 
not as strong as compared with those between Fx and 
Fq, and Fq and Fk (Cheng et al. 2000). Our present 
result is consistent with our previous ones. Besides, we 
wish to mention that there are 10 sources without any 
available redshift. There is no difference between the 
results obtained with those 10 sources both excluded and 
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included. 

Physically, if high-energy electrons in the jet are re- 
sponsible for the radio and 7-ray emission, it would then 
be possible that radio photons are produced by the syn- 
chrotron radiation of these electrons, and that 7-rays are 
produced by the inverse Compton scattering of ambi- 
ent lower energy photons, which may be the synchrotron 
photons. Therefore, there should exist a correlation be- 
tween the radio and 7-ray luminosity to some extent. In 
fact, it is one of the predictions in the synchrotron self- 
Compton (SSC) models (e.g. Maraschi et al. 1992). Ob- 
viously, the seed photons may also be radiated directly 
from the accretion disk or disk photons scattered by a 
broad line region or wind material above the disk into 
the jet path (e.g. Dermer, Schlickeiser 1993; Sikora et al. 
1994; Blandford, Levinson 1995; Zhang, Cheng 1997). 
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